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ABSTRACT 

The recent Kepler discovery of KOI-152 reveals a system of three hot super-Earth candidates that are 
in or near a 4:2:1 mean motion resonance. It is unlikely that they formed in situ, the planets probably 
underwent orbital migration during the formation and evolution process. The small semimajor axes of 
the three planets suggest that migration stopped at the inner edge of the primordial gas disk. In this 
paper we focus on the influence of migration halting mechanisms, including migration "dead zones", 
and inner truncation by the stellar magnetic field. We show that the stellar accretion rate, stellar 
magnetic field and the speed of migration in the proto-planetary disk are the main factors affecting 
the final configuration of KOI-152. Our simulations suggest that three planets may be around a star 
with low star accretion rate or with high magnetic field. On the other hand, slow type I migration, 
which decreases to one tenth of the linear analysis results, favors forming the configuration of KOI- 
152. Under such formation scenario, the planets in the system are not massive enough to open gaps 
in the gas disk. The upper limit of the planetary masses are estimated to be about 15, 19, and 
24 M®, respectively. Our results are also indicative of the near Laplacian configurations that are 
quite common in planetary systems. 

Subject headings: (stars: KIC 8394721) planetary systems: formation-solar system: formation-stars: 
individual (KOI-152) 



1. INTRODUCTION 

The Kepler Mission, launched in March 2009, photo- 
metrically monitors a large patch of sky with sufficient 
precision to detect terrestrial sized planets in potentially 
habitable orbits. The Kepler Mission has released their 
first 16 months data of 2321 transiting pla net candidates 
(jBatalha et al.ll2012tlFa"brvckv et al.ll2012j ). The mission 
is sensitive to a larger rang e of semimajor axe s than 
ground-based transit surveys (|Borucki et all2010ft . Thus 
there are opportunities to detect multiplanetary systems. 
According to statistical results on the first four months 
data, ~ 17% system s have multiple planet candidates 
l|Borucki et al .1120 111 ). More than a dozen multiple tran- 
siting planet systems will be confirmed or rejected by 
mean s of transit timing variations (TTVs) (jFord et al.1 
1201 If ). The candidate multiple planetary systems show 
that at least ~ 16% contain a pair of planets close to 
2:1 period commensurability (with a period ratio of tw o 
planets ranging from 1.83 to 2.18) (|Lissauer et al.ll2011[ ). 
Furthermore, the detailed statistical analysis of the re- 
lease data over 16 months suggests that the fraction of 
multiple planetary systems has raised from 17% to 20% 
and the number of the resonant systems also increases 
(|Batalha et al.ll2bl2l ). The resonant systems are of value 
to the researchers, who investi gate formation and e volu- 
tion of the planetary systems. ISteffen et al.l (|2010f ) ana- 
lyzed five Kepler target stars and their planet candidates, 
in which KOI-152 (KOI: Kepler Objects of Interest) is a 
known system, consisting of three planetary candidates. 
Moreover, three planets, orbiting about an F dwarf star 
(KIC 8394721) with a mass of 1.4 M©, are very close 
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to 4: 2:1 mean motion resonance (MMR) (jSteffen et al.l 
l2010f ). Table 1 shows the orbital elements of the planets 
in some detail, i.e., 152.01, 152.02, and 152.03 in the or- 
der that the transit detection software identified them in 
the Kepler data, and hereafter we label them as Planet 
01, 02, and 03, respectively, for brevity. Their masses are 
now estimated to be in the range of (20-100) M®, (9-30) 
M®, and (9-30) M®, respectively. Considering the den- 
sity limit and formation scenarios, observers may assume 
them to be 60 M®, 15 M®, and 15 M®, respectively. Al- 
though the semimajor axis of Planet 01 is not well deter- 
mined, estimated value indicates that they are very close 
to the central star. The orbital period ratios of each 
pairs are P02/P03 ^ 2.033 and P 01 /P 02 > 1.896. The 
eccentricities are estimated to be zero an d no evidence 
for la rge eccentricities has been revealed (ISteffen et al.1 
[2010h . Therefore, we assume the eccentricities to be zero 
in this work. Based on the possibility of observing such 
a three-planet system, the planets likely occupy nearly 
coplanar orbits with a sm all deviation of incl ination from 
fundamental framework (|Steffen et al.|[2010l ). 

It is well-known that the planetary configurations in- 
volved in 4:2:1 MMR, also occur both in our solar system 
and exoplanetary systems. For instance, the Galilean 
moons of Jupiter had been revealed in a three-body 
Laplace resonance over several hundred years. In ad- 
dition, another example is that three super-Earths con- 
stitute the HD 4 0307 system, where three pla nets are 
near 4:2:1 MMRs (jPapaloizou fc Terquemll2010f ). similar 
to KOI-152. For HD 40307, three proto-planets formed in 
the proto-planetary disk with configuration near MMRs, 
which strongly constrains the planetary formation and 
orbital migration theories. In this sense, the formation 
scenario of such near Laplacian configuration is of great 
interest to the researchers, and such investigations may 
test some planetary formation theory. 

According to core-accretion model, a planet with semi- 
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major axis a ca n grow up with t he material surrounding 
it to mass mi so (jlda fe Linll2004l ) 

M niR ,3/2 3/2/ Afc \ 3/2 / a \ 3/ V M *V 1/2 Af 
M iso * 0.16// 7i / c ^_ j (— ) ^— ) M e 

(1) 

where = (m/3M*) 1 / 3 a is the Hill radius of a planet 
with a mass m, and Af z ~ 7 — 10 i?h is the so-called 
feeding zone of the planet, is the mass of central 
star, /d is the enhancement factor to the Minimum Mass 
Solar Nebula (MMSN), 7; cc is the volatile enhancement 
for exterior or interior to the snow line aj CC with a value of 
4.2 o r 1, respectively. For KOI-152, we adopt ai C0 = 5.29 
AU (|Ida fc Lin! I2004D . If three planets formed in the 
region as shown in Table [T] with their estimated minimum 
masses, fd should be at least 37. Hence, we may conclude 
that it is impossible that all planets formed in situ. 

Now, it is widely believed that there are mainly two 
formation mechanisms to produce short-period planets, 
e.g., planet-planet s cattering and planetary migration 
(|Rasio fc Fordlfl996[) . The planet-planet scattering sce- 
nario always requires a massive planet to stir up the 
eccentricities of other bodies and trigger the scattering 
process. After this proc ess had done, some planets may 
obtain high eccentricity ([Rasio away from 

their initial locations. However, the planets of KOI-152 
are near 4:2:1 MMR with nearly circular orbits, and in 
this sense it seems to be impossible to tune scattering 
scenario to yield such a precise configuration. 

Another mechanism is t hat the orbital migration oc- 
curs in the gas eous disk (iGoldreich fe Tremaind [19801 : 
iLin et all I1996D . Multiple planets are likely to be 
captur ed into MMRs given an appropriate migration 
speed (jMasset fe Snellgrovel[200l . As MMR commonly 
emerges in the planetary systems, the orbital migration 
is now considered as one of the plausible mechanisms to 
form such systems. Take the GJ 876 system as an exam- 
ple, consisting of four planets, GJ 876 b and GJ 876 c 
are locked into 2: 1 MMR, which can be explained by mi- 
gration scenario (iLee fe Peale||2002l; Ui et al.ll2002l 120031: 
IZhou et al.H2005t Izhang et al.lhoiof ). 

Recently, an alternative scenario of collision-merger 
has been propo sed to account for the formation of short- 
period planets (jJi et al.ll2011l ). In this mechanism, sev- 
eral embryos can be excited by giant planets after the 
gas of the disk depletes, and then merge into a larger 
body moving on a close-in orbit. Nevertheless, in these 
scenarios, planets are not easy to form in MMRs. 

Considering the aforementioned factors, we propose a 
scenario to produce a configuration of KOI-152. First, 
three planets are assumed to have formed in the region 
away from the star with their nominal masses. Then, 
the planets may undergo type I or II migration due to 
interaction with the gaseous disk until they halt at the 
inner region of the disk. In this phase, three planets are 
trapped into MMRs during the migration. Finally, tidal 
effects, arising from the central star, circularizes their 
orbits. 

Such a scenario had been supposed by 
iTerquem fe Papaloizoul (|2007T ) . They calculated a 
series of runs that are composed of 10-25 planets or 
planetary cores in a disk with masses ranging from 
0.1 to 1 Mg, which undergo type I migration. They 
showed that hot super-Earths or Neptunes do not 



become isolated during their inward migration, and 
the companions on near-commensurate orbits always 
survive. Nevertheless, in their work, the authors did 
not consider the effect of different speed of type I 
migration, which can o ccur. In a more recent work, 
iPapaloizou fe Terqueml (|2010( ) investigated the forma- 
tion of the HD 40307 system, bearing a resemblance to 
KOI-152. They adopted a similar formation scenario 
and found that in the end of the simulations, the planets 
are driven out of Laplace resonances due to tidal effects 
with the central star, and they finally reach a planetary 
configuration v ery c lose to HD 40307. Furthermore, 
IWang fc Zhoul (|2011l ) accounted for the contribution of 
type I migration and mainly focused on the speed of 
migration under the perturbation of gas-giants in the 
outer region in M dwarf system. They further point out 
that pairs of planets favor forming not only near the 2:1 
MMR but other first order resonances. 

One of the important factors that influence final con- 
figurations is holes in the gas disk. In general, holes 
can appear at two positions. One is the boundary of 
the dead zone and the active zone in the midplane of 
the disk responding to the magnetorotational instability 
MRI) . For protostellar gas disk in ad hoc a-prescription 
Shakura fc Sunvaev I 1 19731) . the mass accretion rate, 
M g = 2>irac s hYi g , is constant across the whole disk, where 
h = c s /Qk is the disk scale height, c s , fi^, and S s are 
the sound speed of the midplane, the Kepler angular ve- 
locity, and the gas density, respectively. As the value of 
a decreases from the active zone (~ 0.018; Sato et al. 
2000) to the dead zone (~ 0.006), in the midplane the 
value of the gas density in the dead zone (outer) is three 
times of that in the active zone (inner). Then a maxi- 
mum density location occurs. Another is the inner hole 
of the disk caused by the coupling of the star's magnetic 
field with the gas. Gas falls onto the surface of the central 
star under the effect of the torque induced by the stellar 
magnetic field, a truncation happens at the inner region 
and an inner hole appears. A maximum density of the 
gas disk appears at the boundary of the inner hole. Due 
to the variation of density, the speed of type I migration 
may be changed. The two mentioned regions play a sig- 
nificant role in forming the final planetary configuration, 
especially for low-mass planetary systems. 

In this work, we focus on exploring the configuration 
formation of KOI-152, mainly on the following aspects: 
(a) the speed of type I migration of planets, (b) the den- 
sity profile of gas disk, (c) the possible range of masses 
of three planets, and (d) the nature of the star in the 
system. This paper is organized as follows: In §2, we 
introduce the adopted disk model, orbital migration and 
eccentricity damping models in the investigation. In §3., 
we present the simulation method and outcomes. Finally, 
we conclude and summarize the results in §4. 



2. MODELS 

2.1. Disk Model 

In order to explore the configuration formation, we 
consider a system consisting of a central star and three 
planets, which formed far away from the star. We assume 
that three planets are initially embedded in a gaseous 
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TABLE 1 

Orbital parameters of the KOI-152 planetary system. 



ID 


Scmimajor axis 


Period 


Eccentricity 


Mass 




(AU) 


(day) 




(A%) 


152.03 


0.124 


13.48 


0.00 


15 


152.02 


0.199 


27.40 


0.00 


15 
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0.305 


> 51.94 
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Fig. 1. — The density profile of the gas disk and the values 
of (5. The upper panel (a) shows the gas density profile with star 
accretion rate (M = [1 x 10~ 9 , 5 x 10~ 8 ] M Q yr _1 ), a doad = 0.001, 
and oj mr i = 0.01. The parameters of the star are set to be B* =0.5 
KG and ii* = 2.5 Rq except the red line for the cases in Group 
4 with high magnetic field B* =2.5 KG. The bottom panel (b) 
illustrates the values of /3, taking M = 2 X 10~ 8 Mq yr — 1 as an 
example. Two maximum density locates at ~ 3.6 days and 26 
days, respectively, where changes from positive to negative at 
DM2 and the region inside DM1. 
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Fig. 2. — The relationship between the gas density and star 
accretion rate. The upper panel (a) shows the orbital period at 
the location of maximum density with various star accretion rate. 
The black dot line represents the boundary of the dead zone and the 
active zone using equation {5{ while the purple line indicates the 
truncated location according to equation I0. The bottom panel (b) 
displays the values of density at two maximum density locations. 



disk. The surface density is given as (|Pringlelll98l 



M 



37ri/(a) ' 



(2) 



where M is the accretion rate of the star and v(a) is 
the effective viscosity at the orbit of a semimajor axis a. 
According to the observation of young cluste r p-Oph, the 
accretion rate of the star can be written as (|Natta et al.l 
L'lHKi: iVorobvov fc Basul[2009l) 



M ~ 2.5 x 10" 



1.3±0.3 



M Q yr- 1 . 



(3) 



According to equation (3), the accretion rate of star 
is ~ 3.87 x 10~ 8 M Q yr -1 for this system. Nevertheless, 
the value will decrease, on average, with the evolution 
of T Tauri star and its disk. Hence, herein we con- 
sider the star accretion rate in the range of [1 x 10~ 9 , 
5 x 10 -8 ] Mq yr -1 . The effective viscosity of the disk is 
v(a) = ac s h, where a and c s , represent the efficiency fac- 
tor of angular momentum transport, and sound speed at 
the midplane, respectively; h = c s /VL means the isother- 
mal dens ity scale height, VI refers to the Kepler angular 
velocity (jShakura fc Sunvaev 1119731 ). Because of the ef- 
fect of MRI, the values of a in the active zone and the 
dead zone are quite different. The effective value of a is 



expressed as (jKretke fc Linll2007t iKretke et aI1 l2009) 



a e g(a) 



erf 



— Qcrit 



0.1a, 



crit 



( 4 ) 

where a mr i , adead denote the value of a in active zone and 
dead zone, respec tively. Herein, we choose adcad = 0.001 
and a mrl = 0.01 (jSano et al.ll2000t ). The parameter a crit 
in the error function erf is the location of the boundary 
of MRI and 0.1a cr it repres ents the width of th e transition 
zone. a cr it is modeled as (jKretke et al.ll2009l ) 

/ ■ \ 4 / 9 



VO02/ 



Mq yr- J 

1/5 



1/3 



lcm 2 g 



(5) 



where kd is the grain opacity. 

Considering the disk depletion, the gas density profile 
can be modified to be 



9 ?>TTv{a) 



exp 



t 

^dep 



(6) 



where T^ ep refers to the disk de pletion timescale, w hich 
is observed as few million years (|Haisch et alll200lD . We 
adopt Tdop = 10 6 yr in the simulations, where t is the 
time of evolution. 
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Because of the s tellar magneti c field, the gas disk is 
truncated at a ms t r (jKocnigl 1991J) 



12/7 



"""" (1JH) 10 x{ '^J (loMc) 



4/7 



-1/7 



M 



-2/7 



10" 7 M Q yr- 1 



(7) 



where i?©, and £?* refer to the radius of the star, the 
radius of the sun and the magnetic field of the central 
star, respectively. j3' < 1 is a free parameter. Herein, we 
choose P' — 1 corresponding to a typical Alfven radius 
in the way of spherical accretion. Hereafter we use P C rit 
and PmstD the Keplerian orbital period, instead of a cr i t 
and a ms tr for convenience. Combining the effect of the 
magnetic field, the gas density profile is substituted by 



M 



E f = q — T\ ex P I I *!> 



t 



where 



r} = 0.5 



^dcp 



. CI Q-mstr \ , ., 
erf ( + 1 

0.1 a ms tr / 



(8) 



0) 



is induced by the truncation of the magnetic field. 

Based on equation ©, we find that S s oc r _1 and 
there are generally two maxima in the gas density pro- 
file. Figure Q] shows the density profile with various star 
accretion rates (see the top panel (a)). We label the inner 
maximum density location as DM1, similarly, the outer 
one as DM2. The locations of DM1 and DM2 change 
with the star accretion rate. In Figure [2] we show the 
density versus star accretion rate. From Figure [2] (the 
top panel (a)), we notice that, with a decrease of star ac- 
cretion rate, the values of DM1 and DM2 approach each 
other until they merge into one, the combination occurs 
at ~ M = 5 x 10~ 9 M© yr -1 , corresponding to a maxi- 
mum density at an orbital period of 13 days. Herein we 
choose ctdoad = 0.001, a mi i = 0.01, and £>* = 0.5 KG, 
respectively, where the bottom panel (b) displays the val- 
ues of the density at DM1 and DM2, respectively. 

2.2. Planetary migration and eccentricity damping 

For the planets in KOI-152 system, their masses are 
estimated to be less than 100 M®. From a classical plan- 
etary formation theory, they may undergo type I or type 
II migration during their evolution. 

Type I migration is induced by angular-momentum 
exchange between gas disk and planets. Based on lin- 
ear analysis, the net momentum loss on a planet causes 
an inward migration dGoldreich fc Tremainelll979t IWardl 
Il997t iTanaka et "all l2002t) . Under this assumption, the 
speed of type I migration is very fast. In such situa- 
tion, it is difficult to produce terrestrial planets. Re- 
cently, there are several new theories on reducing the 
speed of type I migration or even reverse the migrat- 
ing direction dBaruteau fe Massetl 120081: Kiev fe Cridal 
2008tlKlev et al.l2009HPaardekooper fe Papaloizoul2008t 
Wang fe Zhoull201lD . 

Considering the uncertainty of type I migration, we 
adopt a reduction of the migration speed, taking a 
timescale of type I migration of an embryo with mass 



to as (jTanaka et alJl2002T > 
1 



a 



Tmkrl 



-TIL 



1 (M* 
/i(2.7+ 1.1/3) 





"1 4. ( er \5- 




1 ( er \4 
I 1 - \l.lh> J 



rrV, (io) 



where e, r, h, and f2 are eccentricity, distance from the 
central star, scale hight of the disk, and the Keplerian 
angular velocity, respectively. TLinear is the timescale of 
linear analysis result, /i is the reduction factor. Herein, 
we choose fi — 0.03, 0.1, and 0.3, respectively, and 



-dlnE g /dlna. T, g means the gas density profile 



of the disk expressed in equation ((HJ. As the value of /? 
is related to the density profile, the speed of type I migra- 
tion may be slowed down or even reversed at some special 
areas. In addition, Figure 1 shows the values of j3 using 
M = 2 x 10~ 8 Mq yr -1 as an example (see the bottom 
panel (b)). According to equation (10), if < —2.45, the 
timescale of type I migration will transfer from positive 
to negative. Furthermore, we notice that when embryos 
run across DM2 and DM1, the migration speed will de- 
crease, which may lead to the trapping of the embryos 
there. 

When it grows massive enough, a planet will start 
to experience type II migration, as the strong torque 
caused by the planet wil l open a gap in the gaseous disk 
(|Lin fc Papaloizoulll993D . The time scale of type II mi- 
gration for a planet with mass to is (|Ida fc Lidl2008h 

TmigUl ~ 5 X lO 5 /^ 1 

( , , , , 

VlAU 



Vio- 



M, 



1/2 



s-112 



-7 = 0.7 x 10- 

a 



V1O- 3 / VlAU/ M fl 



Mq 

M, 



-1/2 



yr, 



-1/2 



yr,( 



where a is the efficiency factor of angular momentum 
transport. When the mass of planet is comparable to 
that of the gas disk, a fraction (C2 ~ 0.1) of total angular 
momentum will transfer between the planet and the disk 
in the evolution. In this case, r m i g ni applies. Herein 
we adopt a = ayead = 0.001, Cict = 10~ 4 . If the mass 
of the gas disk is larger than that of the planet, it will 
migrate with the gas disk over the timescale of r m i g n2- 
We emphasize that the timescale of type II migration is 
larger than that of type I. 

A gap will form in the gas dis k when the plane t grows 
to massive enough (to > M crit ) (jlda fc Linll2008| ). 

— °(i^)(rlu) 1/2 (^) M - (I2) 

In this work, we assume that the planet will undergo type 
II migration when its mass is greater than M cr it- 

Additionally, we consider the eccentricity damping in- 
duced b y the interactions between the gas disk and 
embryo (|Goldreich fc Tremaind Il979f ). The damping 
timescale for an embryo with mass m is described as 
(|Cresswell fc Nelsonll2006h 



(I) 

\ 6 / edamp 



Qe / M* 



0.78 V m 



1 / r 
4 V h 



2^ r 



(13) 
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where Q e = 0.1 is a normalization factor to be consis- 
tent with hydrodynamical simulations. The meanings of 
other symbols are similar to those in equation (fTU|) . 

3. NUMERICAL SIMULATIONS AND RESULTS 

To explore the secular evolution of the KOI-152 plan- 
etary system, we assume that three planets had formed 
in the outer region of the system and migrate toward the 
inner region due to the interactions with the gas disk. 
Thus, the acceleration of a planet with mass rm is given 
as 



G(M» +m 4 ) 



N 

£ 



where and mean the position and velocity vectors 
of the planet vrti in the stellar-centric coordinates, and 
the external forces are defined as 



cdamp 



migl 



migll 



,(Yi 



rfr e 



^Tmigl ' 



(15) 



2T m igH 

Each planet is assumed to be initially in coplanar and 
near-circular orbit and suffers from mutual gravitational 
interaction and the effect exerted by the central star. The 
initial orbital elements of each planet were randomly gen- 
erated: argument of pericentre, longitude of the ascend- 
ing node, and mean anomaly were randomly set between 
0° to 360°. In this way, we generate a set of 34 runs for 
simulation. 

We numerically integrate the equations (1141 using a 
time-symmetric Hermit scheme (jAarsethll2003h . If the 
distance from a planet to the central star is smaller than 
the radii of the star we originally set, we assume that the 
planet collides with the star. Each run evolved for 10 s 
yr. 

Based on equation ([P2]). for a = 10~ 3 , if a > 0.125 
AU, then we have M cr j t > 15 M®. Hence, if their ini- 
tial locations are all without 0.125 AU, KOI-152. 02 and 
KOI-152. 03 will undergo type I migration, the outmost 
one will experience type I or II migration. In this sense, 
two kinds of models are assumed in the simulations: For 
Model 1, we suppose that all planets will suffer type 
I migration; whereas for Model 2, we simply consider 
that KOI-152.02 and KOI-152.03 will go through type 
I migration but KOI-152. 01 will undergo type II migra- 
tion. According to the aforementioned analysis, one of 
the planets in this system may be captured at DM1 or 
DM2. In such case, taking into account the range of the 
star accretion rate [1 x 10~ 9 , 5 x 10~ 8 ] M yr -1 , we find 
that the positions of DM2 are at ~ 13, 26, and 48 days, 
respectively, in three groups for Model 1. In addition, 
the conditions of density profile of various groups and 
models are summarized as follows. 

Model 1 

Group 1: For M = 5 X 10~ 8 M & yr^ 1 , the density 
profile S 5 is labeled as the blue line in Figure [T] DM2 = 
48.039 days, KOI-152. 01 is likely to be captured at a cr it- 



1 1 r 


1 1 1 1 

152.03 

152.02 

152.01 


i ! 

(a) i 


I- 












1 


i . i 





F mig i(or + F mig ii), (14) 



(b) 



1.0 

Time (Myr) 



Fig. 3. — The evolution of orbital periods and resonant angle in 
Group 1. At ~ 0.2 Myr, KOI-152.02 and KOI-152.03 are evolved 
into a 2:1 MMR in the evolution. The upper panel (a) shows the 
time evolution of the orbital periods. Planet 03 halts at DM2 
first until Planet 02 comes. Planet 02 kicks 03 into the inner re- 
gion when it migrates towards DM2 and the outmost planet keeps 
Planet 02 repeating this process just as the first two planets do. At 
the end of the run, Planet 03 stops at DM1. The bottom panel (b) 
shows one of the 2:1 resonant angles of Planet 02 and 03. (Here- 
after, we label 8\ = 2Ao2 — A03 — ^02, #2 = 2A02 — A03 — roo3i 83 = 
2Aoi — A02 — rooi, 84 = 2Aqi — A02 — roo2-) 



Group 2: For M = 2 x 10~ 8 M Q yr" 1 , E g is labeled as 
the black line (FigureE]). DM2 = 26.08 days, KOI-152.02 
is likely to be captured at a cr it. 

Group 3: For M=lx 10" 9 M Q yr" 1 , E ff is labeled 
as the green line (Figure [J). DM2 = 13.067 days, KOI- 
152.03 is likely to be captured at a cr it- 

Group 4: For M = 1 X 10~ 9 M Q yr" 1 , E g is labeled as 
the red line (Figure HJ). DM2 = 51.9 days, KOI-152.01 is 
likely to be captured at a cr j t . 

Model 2 

Comparison with Group 2 of Model 1: For M — 
2 x 10~ 8 Mq yr -1 , E g labeled as the black line in Figure 
[TJ DM2 = 26.08 days, KOI-152.02 is likely to be cap- 
tured at cj cr it. However, KOI-152.01 will undergo type II 
migration. 

3.1. Model 1: All planets undergo type I migration 

In order to understand the influence of the star ac- 
cretion rate and the speed of type I migration, we per- 
form four groups of simulations. Table 2 are listed the 
detailed information for the dominatant results of each 
group, where the last column shows the Keplerian orbital 
periods at the locations of DM1 and DM2. 

3.1.1. Group 1: KOI-152.01 captured at a cri t 

In this group, DM1 and DM2 occur at <~ 2.4 days and 
48 days, separately, with M = 5 x 10~ 8 M Q yr -1 . From 
the profile of gas density, we note that, if Planet 01 is 
likely to be trapped at DM2, the resulting configuration 
may be analogous to KOI-152. In total, we performed 
five runs to examine how the planets come into the res- 
onant region, where Table 2 reports typical outcomes in 
our simulations. For example, in a typical run, Planet 02 
or 03 will continue to move inward until reaching DM1 in 
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Fig. 4. — The results of the evolution for Group 3-2. The two upper panels (a and b) show the evolution of orbital periods and 
eccentricities. The middle and bottom panels (c-f) are indicative of the resonant angles of three planets. In this run, they are captured 
into 4:2:1 MMR over very shorter timescale. In the end of this run, a planetary configuration similar to KOI-152, which consists of three 
planets, is formed. 



TABLE 2 

Detailed information of five groups in the simulations. Terminal periods are the orbital periods of three planets at the 
end of run. p cr ]t and p ms tr represent keplerian orbital periods in association with the regime of dm1 and dm2. 



ID 


M 


Initial periods 


fl 


Terminal periods 


Terminal periods ratios 






(M© yr" 1 ) 


(days) 




(days) 


(P01/P02, P02/P03) 


(days) 


Group 1 


5 x 10" B 


120, 320, 850 


0.1 


2.58, 5.2, 60.5 


11.63, 2.02 


48.039, 2.443 


Group 2-1 


2 x 1CT 8 


120, 320, 850 


0.03 


3.8, 7.75, 32.83 


4.24, 2.04 


26.08, 3.618 


Group 2-2 


2 x icr 8 


120, 320, 850 


0.1 


1.899, 12.19, 24 


1.97, 6.42 


26.08, 3.618 


Group 2-3 


2 x icr 8 


120, 320, 850 


0.3 


3.8, 7.67, 32.85 


4.28, 2.02 


26.08, 3.618 


Group 2-4 


2 x icr 8 


120, 320, 20012 


0.03 


31.72, 64.23, 6960.77 


108.37,2.02 


26.08, 3.618 


Group 2-5 


2 x icr 8 


120, 320, 20012 


0.1 


3.82, 7.7, 32.85 


4.27, 2.01 


26.08, 3.618 


Group 2-6 


2 x icr 8 


120, 320, 20012 


0.3 


3.78, 5.72, 32.83 


5.74, 1.51 


26.08, 3.618 


Group 3-1 


l x icr 9 


120, 320, 850 


0.1 


3.5, 7.12, 14.73 


2.07, 2.03 


13.067 


Group 3-2 


l x icr 9 


70, 150, 320 


0.03 


13.78, 27.87, 57.07 


2.05, 2.02 


13.067 


Group 3-3 


l x icr 9 


220, 500, 1400 


0.1 


11.71, 23.56, 47.64 


2.02, 2.01 


13.067 


Group 4 


1 x 10" 9 


220, 320, 850 


0.1 


14.17, 28.56, 58.54 


2.05, 2.01 


51.9 


KOI-152 








13.48, 27.40, 52.09 


1.90, 2.03 





cases where they have been kicked inside of DM2. Sub- 
sequently, Planet 01 is trapped at ~ 60 days as shown in 
Figure 3 and the two inner planets are eventually cap- 
tured into 2:1 MMR. From Figure 1, we show that if a 
planet is pushed into the region inside DM2, the value of 
f3 changes from negative to positive before it approaches 
DM1. The change in sign of /3 is the reason that the two 
inner bodies migrate towards DM1. However, no matter 
how we varied the speed of type I migration or re-scaled 



the initial position of Planet 01, (e.g., much more distant 
than the value given in Table 2), Planet 02 and 03 still 
rush into the inner zone about DM1 in the evolution. 
Hence, we may safely conclude that a configuration like 
KOI-152 cannot be generated when the star is too young 
with high star accretion rate. To sum up the above out- 
comes, for five runs, we find that Planet 02 and 03 are 
locked into 2:1 MMR in three cases during dynamical 
evolution, whereas they are captured into 3:2 MMR in 
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Fig. 5. — The results of the evolution Group 3-3. The two upper panels (a and b) show the evolution of orbital periods and eccentricities. 
The middle and bottom panels (c-f) exhibit their resonant angles. In this run, they are captured into 4:2:1 MMR within a shorter timescale. 
A planetary configuration analogous to KOI-152, consisting of three planets, is generated. 



1600 | r 




350 



Orbital Period (day) 

Fig. 6. — The gas density profile in Group 3 and 4. Because of 
a high magnetic field of the star, the gas disk will be truncated at 
distant region with low gas density in Group 4. 



the other two runs. 

In this scenario, a planetary configuration is finally 
created, that consists of two planets trapped into a 2:1 
MMR or 3:2 MMR at DM1 and the outmost one resided 
inside DM2, far away from DM1. 



3.1.2. Group 2: KOI-152.02 captured at a cri t 

By adopting M = 2 x 10~ 8 M yr _1 , we have DM1 
and DM2 located at ~ 3.6, 26 days, respectively. In 
this group, we performed six runs in total. The results 
of all runs are reported in Table [2] From Table [2j we 
observe that, similar to the cases of Group 1, Planet 03 
is quickly captured at DM1 as it is thrown into the re- 
gion inside DM2. In Run 1-3, we utilize the same initial 
conditions but choose a variational migration speed. In 
these runs, Planet 02 cannot stay at DM2 but is kicked 
into the inner region by the perturbation of Planet 01, 
then it falls into MMR with Planet 01 (run 2) or Planet 
03 (run 1,3). In Run 4-6, Planet 01 is initially well sep- 
arated from the other two, by placing the outer planet in 
a starting location of much more distant from the central 
star. The results of Run 5 and 6 are analogous to that 
of Run 1-3. The variation of /3 is the governing reason 
for the evolution of Run 1-3 and 5-6. However, owing 
to slow migration speed, Run 4 differs from other runs, 
where three planets cannot pass through DM2. Here, 
one may notice that it is impossible for three planets to 
form a configuration resembling to KOI-152. From the 
analysis of six runs, we summarize the simulation out- 
comes, where Planet 02 and 03 are in a 2:1 MMR for 
four runs and in one run they are captured into a 3:2 
MMR , whereas in another special case Planet 01 and 02 
become trapped into a 2:1 MMR over the evolution. 

The created configuration in Group 2 varies with the 
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Fig. 7. — The results of the evolution for Group 4. The magnetic field is B, = 2.5 KG. Note that three planets are captured into 4:2:1 
MMR. The middle and bottom panels (c-f) show the time variations of their resonant angles. 



value of f\. The inner planet is always trapped in the 
boundary of the inner hole unless they undergo a lower 
speed of type I migration and the outmost planet is 
formed farther simultaneously. Thus, in such scenario, 
a configuration similar to KOI-152 cannot be created. 

3.1.3. Group 3: KOI-152. 03 captured at a crit 

In the simulations, we assume the star accretion rate 
is equal to 1 x 10 -9 M Q yr -1 . From the density profile 
shown in Figure [TJ we see that DM1 and DM2 tend to 
combine into one position at ~ 13 days. In this group, we 
totally carried out 10 runs. Details of three runs are also 
shown in Table [2j If we examine the type I migration 
considering fx > 0.1, the results of evolution then bear 
resemblance to the case of Group 3-1, where KOI-152. 01 
is trapped at ~ 14 days in the meantime the other two 
planets jump into the inner region at P < 10 days, unless 
we started them at distant orbits from the star or slow 
down the speed of type I migration. In Group 3-2, we 
lower the migration speed and the result is shown in Fig- 
ure 21 where Planet 03 is capable to being trapped at ~ 
14 days. The simulation results are consistent with the 
current observational values of KOI-152. When Planet 
02 and 01 keep pace with Planet 03, two pairs are all 
captured into 2:1 MMRs. Group 3-3 simulates the con- 
dition related to a high speed of type I migration and 
distant initial orbits. Figure [5] displays that Planet 01 
and 02 are locked into 2:1 MMR in the dynamical evolu- 
tion as well as Planet 02 and 03. In conclusion, according 



to our evaluation, three planets are eventually captured 
into Laplace resonance, with the resonant angles in a li- 
bration about 180° for Group 3-2 and 3-3. When the 
planets are trapped into MMRs, their eccentricities are 
excited during the evolution; but due to strong gas damp- 
ing, they cannot be pumped up to large values. 

Taking j\ < 0.1 into account, a configuration similar 
to KOI-152 is formed, where two pairs of planets are in 
the 2:1 MMRs. 

3.1.4. Group 4 : High magnetic field star 

Herein we consider that the star bears a high magnetic 
field of = 2.5 KG. In comparison with Group 3, for 
a high magnetic field star, the truncation of the inner 
hole is at ~ 50 days, farther than that a low magnetic 
field star as shown in Figure 6. Herein we carried out 
three runs. Figure 7 illustrates the results of a typical 
run, where the middle and bottom panels (c-f) show the 
resonant angles for the planets. From the figure, we note 
that three planets are captured into a 4:2:1 MMR at ~ 
1 Myr. Simultaneously, we examine the resonant angles 
and find that three planets are also locked into a Lapla- 
cian resonance, librating about 0° at ~ 4 Myr. When 
they are captured into MMR, Planet 02 and 03 may mi- 
grate towards the inner region of the hole with little gas. 
Consequently, the influence of the gas damping is not 
dominant. Thus, the eccentricities of Planet 02 and 03 
are stirred up to ~ 0.2 when they are trapped into a 
2:1 MMR. However, Planet 01 remains outside the hole, 
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middle and bottom panels (b and c) exhibit the time evolution of 
the semimajor axes and eccentricities, respectively. 



and its eccentricity finally decreases to ~ 0.002. For one 
paradigm, because of fast speed of type I migration, the 
system is totally destroyed. In another run, the produced 
planetary configuration is quite analogous to the typical 
Laplacian geometry, with the resonant angles librating 
at 180°. For fi < 0.1, we have a planetary configuration 
similar to KOI-152. 

3.2. Model 2: KOI-152. 01 undergoes type II migration 

According to our above analysis, Planet 01 may per- 
form type II migration. In this model, we set Planet 02 
and 03 to undergo type I migration but Planet 01 suffers 
from type II migration. We carry out 10 runs for Model 
2, by varying the initial position of the outmost planet. 
But owing to a slow speed of type II migration, Planet 
01 cannot reach its nominal position even if it is placed 
at the orbit with a period of 70 days in the beginning. 
Figure [8] shows a typical run of the simulations. In this 
run, the density profile is the same with that in model 
1 Group 2, /i=0.1 for Planet 02 and 03. In the end, 
the simulation results show that three planets migrate 
to 31.72, 63.86 and 211.13 days, respectively. However, 
Planet 01 is a little deviated from the present-day obser- 
vation. 

Under this circumstance, Planet 01 cannot approach 
its estimated position regardless the initial location of 
the outmost planet. In addition, the inner two planets 
are always trapped into 2:1 MMR at DM2. In a word, 
we cannot generate a Laplacian configuration for three 
planets from Model 2. 

4. CONCLUSIONS AND DISCUSSIONS 

In this work, we have extensively investigated the for- 
mation of configuration for the KOI-152 system using 
numerical simulations. We assume three planets formed 
in a region far away from their current locations. The 
embedded planets in the gas disk will migrate into the 
regime much closer to the central star. In order to pro- 
duce a configuration similar to KOI-152, there are some 



requirements for the disk and the star. In summary, we 
reach the following conclusions: 

1. If KOI-152. 01 survives as predicted, no type II mi- 
gration happened for three planets in the system. 

2. From equation (12), the masses of three plan- 
ets can be limited to be mo3 < 15, mo2 < 19, and 
moi < 24 Mq, respectively. According to core-accretion 
model, a solid core may continue to ac crete gas en- 
velop until it grows up i nto 2 - 10 M m (iBarnes et al.l 
2009t iBodenheimer et ail 120001: iHubickvi et al.l 120051: 
Ikoma et al.l 120011 : iPollack et al.lfl996|) . therefore simply 
judging from the radii of the planets and the formation 
scenario, we cannot clearly identify the composition of 
the planets, and their masses are not well determined. 
Based on the gap opening criteria, the masses of three 
planets are evaluated to be (9-15) M®, (9-19) M®, and 
(20-24) M®, respectively. 

3. Low star accretion rate and high magnetic field of 
the star are propitious to the formation of the config- 
uration of KOI-152. Low star accretion rate may im- 
ply that the formation of three planets takes place at 
the later stage of the star evolution, where DM1 and 
DM2 are merged into one. From Figure [21 we learn that 
the boundary is about M — 5 x 10 _9 Af Q yr _1 . In this 
sense, one requirement for the formation of KOI-152 is 
M<5X 10- 9 A/ o yr" 1 . 

4. Low speed of type I migration with /i < 0.1 facil- 
itates the formation of the KOI-152 system. This con- 
clusion agrees with those mentioned by Ida & Lin (2008) 
and Wang & Zhou (2011). 

5. In this work, we had not examined tidal inter- 
actions between planets and central star after plane- 
tary formation scenario. On the basis of tidal timescale 
(|Mardling fc Lin|l2004t IZhou fc Linll2008ft . the eccentric- 
ities of Planet 01, 02, 03, with their nominal orbital el- 
ements will be damped in r t id c , 1-3 x 10 9 Q', 2 x 10 8 Q', 
and 9.6 x 10 6 Q' years, respectively, where Q' is the tidal 
dissipation factor and we adopt p = 3 g cm -3 for them. 
Meanwhile, the semimajor axes are also decreasing over 
the timescales of Ttide/(2e 2 ). In that case, the orbits of 
the planets will become a little closer to the star that 
those without consideration of tidal effects. Owing to 
shorter timescales, the inner planets migrate faster, driv- 
ing all of them out of MMR. Then, the configuration 
of KOI-152 forms and three planets are in near MMRs. 
This scenario may be suitable to account for the forma- 
tion of other systems with planetary configurations like 
KOI-152. 

Among 16 months data of Kepler, 242 target stars 
host two planet candidates, 85 with three, 25 with 
four planets, eight w ith five planets and one with six 
(jFabrvckv et all 120121 ). We examine all candidates and 
find that ten systems may have two pairs of the planets 
involved in near 2:1 MMRs, e.g., a near Laplacian reso- 
nance configuration. The orbital periods and their ratios 
are shown in Table 3 . The ratios of orbital periods are 
all in the range of [1.83, 2.18], which may imply that they 
are likely to bear a formation scenario similar to KOI- 
152. In addition, the configuration of KOI-148 is close 
to the outcomes of Group 2-1 and 2-3 where two plan- 
ets are simply trapped in 2:1 MMR, which demonstrates 
that the system may be formed when the star is older 

1 http://archive.stsci.edu 
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TABLE 3 

THE ORBITAL PERIODS OF EACH PLANETARY CANDIDATE AND THE RATIOS BETWEEN THEM. Pol, Po2,-Po3, AND P04 A R E : RESPECTIVELY, THE 

ORBITAL PERIODS OF THE PLANETS IN THE SYSTEMS. 



1JJ 










ri / r.>. ■ 
"01/ "02 


Pi 02/P03 


"03 1 "04 


KOI 


(days) 


(days) 


(days) 


(days) 








152 


52.09119 


27.40415 


13.484 




1.9009 


2.0323 




571 


13.34331 


7.26733 


3.886758 




1.8361 


1.8698 




665 


5.867973 


3.07154 


1.611912 




1.9104 


1.9055 




7.33 


11.34917 


5.924992 


3.132968 




1.9155 


1.8912 




829 


38.5596 


18.64902 


9.75222 




2.0676 


1.9123 




898 


20.08923 


9.77059 


5.16991 




2.0561 


1.8899 




899 


15.36813 


7.11388 


3.306569 




2.1603 


2.1514 




1426 


150.0341 


74.91443 


38.87641 




2.0027 


1.9270 




1860 


12.2094 


6.3194 


3.0765 




1.9321 


2.0541 




1895 


32.1349 


17.2812 


8.4575 




1.8595 


2.0433 




935 


87.6464 


42.6329 


20.85987 


9.6168 


2.0558 


2.0438 


2.1691 


148 


42.89554 


9.67374 


4.777978 




4.4342 


2.0247 





t han the case of KOI -152 with M = 2 x 10 8 M Q yr^ 1 . 

iMarcv et al.l (|2001l) revealed that there were two giant 
planets involved in a 2:1 MMR orbiting the star GJ 876. 
However, after the fourth planet was discovered in the GJ 
876 system, the previously known 2:1 MMR configura- 
tion then becomes a Laplacian resonance configuration, 
where three planets are trapped into a 4:2:1 MMR with 
the resonant angles librating at 0°, respectively, differing 
from the Galilean m oons of Jupiter in a libration at 180° 
([Rivera et al.l l2010). In such configurations, a planetary 
system will remain stable at least one billion years. A 
great many of previous works had investigated the dy- 
namics and stab ility of the sy s tem consisting of three 
planets at then (Ui et al.l I2002L 120031 : IZhou et al.l 120051 : 
Zhan g~et al.l I2010TK From the former results, we learn 
that the migration scenario may be responsible for the 
formation of two giants which is locked into 2:1 MMR. 
Thus under the same formation scenario for KOI-152, 
three planets may be also trapped into a 4:2:1 MMR, 



showing a resemblance to those of GJ 876. Hence, we 
may reach a safe conclusion that the near Laplacian con- 
figurations are quite common in the planetary systems as 
also revealed by Kepler and our work may provide some 
substantial clues to the formation of such intriguing sys- 
tems. 
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